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Abstract 

The vacuum dynamics of SU(2) and SU(3) lattice gauge theories is stud- 
ied by means of a gauge-invariant effective action defined using the lattice 
Schrodinger functional at finite temperature. In the case of the SU(3) gauge 
theory numerical simulations are performed at zero and finite temperature. 
The vacuum is probed using an external constant Abelian chromomagnetic 
field. At zero temperature, in agreement with our previous studies for the 
SU(2) theory, the external field is totally screened in the continuum limit. 
At finite temperature numerical data suggest that confinement is restored by 
increasing the strength of the applied field. We give also an estimate of the 
deconfinement temperature. 
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1. INTRODUCTION 



The lattice approach to gauge theories allows the non-perturbative study of gauge 
systems without loosing the gauge invariance. Recently it has been proposed 
a method to define on the lattice the gauge invariant effective action by using the 
Schrodinger functional 

In the continuum the Euclidean Schrodinger functional in Yang-Mills theories with- 
out matter fields reads: 



1.11 



In Eq. (|1 ■ 1|) H is the pure gauge Yang- Mills Hamiltonian in the fixed-time temporal 
gauge, T is the Euclidean time extension, while V projects onto the physical states. 
^ki^) ^"^^ ^'k^\^) static classical gauge fields, and the state \A) is such that 



1.2) 



being a wavefunctional. Note that, by definition, A*^*^] is invariant 

under gauge transformations of the gauge fields A^^^ and A^-^\ 
Using standard formal manipulations and the gauge invariance of the Schrodinger 
functional it is easy to rewrite ClS db functional integral 



Z[A^^\A^' 



1.3) 



with the constraints: 



A^,{xi = 0) 



A 



1-4) 



v4^(x4 = T)= A 



if) 



It is worthwhile to stress that in Eq. ( |1.3| ) we should sum over the topological 
inequivalent classes. However, it turns out |^ that on the lattice such an average 
is not needed because the functional integral in Eq. ( |1.3| ) is already invariant under 
arbitrary gauge transformation of A^'^^ and A^^\ The lattice implementation of the 
Schrodinger functional Eq. (|1.3| ) is given by [^]: 
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where the integrations over the links U^{x) are done with the fixed boundary con- 
ditions: 



Z[U^,W] = VUe 



U{x)\ 



X4=0 



f/«, U{x)\ 



X4 = 
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The links U^'^ and f/(^) are the lattice version of the continuum gauge fields A^''^ and 
In Ref. in we introduced the new functional: 



where 



and Z[0] means the Schrodinger functional Eq. ( |1.8| ) without external background 
field (f/;;"* = 1). 

From the previous discussion it is evident that ZlW^^^] is invariant against lattice 
gauge transformations of the external link U^^^. Moreover, it can be shown that |l|]: 

lim r[l""*] = Eo[^'"*] - ^o[0] (1.9) 

where i?o[^^^*] is the vacuum energy in presence of the external background field. 
Thus we see that r[A®^*] is the lattice gauge-invariant effective action for the static 
background field A'^^*. It is worthwhile to stress that our lattice effective action 



is defined by means of the lattice Schrodinger functional Eq. (L5) with the same 



boundary fields at 2:4 = and X4 = T. As a consequence we have [Q: 

= / VU e-^"^ , (1.10) 



where Sw is the familiar Wilson action and the functional integral is defined over a 
four-dimensional hypertorus with the "cold- wall" at 0:4 = 0: 

u,{x)u=o = u;-\ (1.11) 

In previous studies |^|-^ we investigated the vacuum dynamics of the SU(2) lattice 
gauge theory by using an external constant Abelian chromagnetic field. In that case 
the relevant quantity is the density of effective action: 



;i.i2) 



where Q = V ■ T, V being the spatial volume. 

The aim of the present paper is twofold. First, we extend the study of the effective 
action with an external Abelian chromomagnetic field to the lattice SU(3) gauge 
theory at zero temperature. Second, we consider both SU(2) and SU(3) gauge 
systems at finite temperature. 

At finite temperature the relevant quantity is the partition function : 

Tr[e-^^^]= j VA{A\e-'^^^V\A) , (1.13) 
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where [3t is the inverse of the physical temperature. The thermal partition function 
Eq. can be written as a functional integral [Q: 



Ja^,(^t,x)=A^,(0,x) 



The lattice implementation of Eq. (|1.14|) is straightforward. We have: 



TT[e-'^^^]= VUf,{xi,x)e-^'^ . (1.15) 

Ju„ (Bt .x)=U„ (0.x)=U„ (x) 



'U^{l3T,x)=U^{0,x)=U^ix) 

Note that, by comparing Eq. ( |1.15|) with Eqs. (|1.10|) , we have: 



Tr [e-^-^] = j VU,{x) Z[U,{x)] , (1.16) 



where Z[Uf^{x)] is the Schrodinger functional Eq. (|1.10|) defined on a lattice with 
L4 = (3t, and "external" links t/^(x) at X4 = 0. 

We are interested in the thermal partition function in presence of a given static 
background field A^^^{x). In the continuum this can be obtained by splitting the 
gauge field into the background field A'^^^{x) and the fluctuating fields ff{x). So that 
we could write formally: 

Zt[A''''] = j Vf] f]\e-'^^^P\ A'''\ff) . (1.17) 

The lattice implementation of Eq. ( 1.17 ) can be obtained from Eq. ( 1.15|) if we write 



Uk{PT,x) = UkiO,x) = f/f* (f)f/fc(f) , (1.18) 

where f/^^*(x) is the lattice version of the external continuum field A'^^^ and the 
Uk{xys are the fluctuating links. Thus we get: 

Zt[A'^'] = I VUkix) VU^ix) Z[U^^\x), Uk{x)] . (1.19) 

Jx4=0 

Note that in Eq. (|1.19|) only the spatial links belonging to the hyperplane 0:4 = are 
written as the product of the external link U^^^{x) and the fluctuating links Uk{x). 
The temporal links [/^{xi = 0,x) are left freely fluctuating. It follows that the 
temporal links (/^{x) satisfy the usual periodic boundary conditions. We stress that 
the periodic boundary conditions in the temporal direction are crucial to retain the 
physical interpretation that the functional is a thermal partition function. 

In the following the spatial links belonging to the time-slice X4 = will be called 
"frozen links" , while the remainder will be the "dynamical links" . From the physical 
point of view we are considering the gauge system at finite temperature in interaction 
with a fixed external background field. As a consequence, in the Wilson action Sw 
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we keep only the plaquettes built up with the dynamical links or with dynamical 
and frozen links. With these limitations it is easy to see that in Eq. ( |1.19| ) we have: 



Ukix) 



z [[/r^x 



;i.2o) 



Indeed, let us consider an arbitrary frozen link f/^^*(x)t/j 
the modified Wilson action by means of the plaquette: 



This link enters in 



PkA{xi = 0,f) = Tr <^ f/f*(f)f/f*(x)f/4(0,f+ A;)f/2(l,f + A;)f/](0,f) \ . (1.21) 



Now we observe that the link t/4(0, x + k) in Eq. ( p. .21] ) is dynamical, i.e. we are 
integrating over it. So that, by using the invariance of the Haar measure we obtain: 



Pm{xa = 0,f) 



Tr|f/i 



cxt 

k 



'x)Ui{Q,x + k)Ul{l, X + A;)f/i(0, x] 



(1.22) 



It is evident that Eq. ( |1.22|) in turns implies Eq. (|1.20|) . Then, we see that in 
Eq. (|1.19|) the integration over the fluctuating links Uk{x) gives an irrelevant multi- 
plicative constant. So that we obtain: 



T 



A 



cxt 



VUe 



-Sw 



;i.23) 



where the integrations are over the dynamical links with periodic boundary con- 
ditions in the time direction. As concern the boundary conditions at the spatial 
boundaries, we keep the fixed boundary conditions Uk^x.Xi) = U^^^{x) used in the 
Schrodinger functional Eq. ( |l.l(j| ). We stress that, if we send the physical tempera- 
ture to zero, then the thermal functional Eq. ( |1.23| ) reduces to the zero-temperature 



Schrodinger functional Eq. (|1.5[) with the constraints Uk{x)\x^=o = Ul^^{x) instead 



of Eq. (|1.11| ). In our previous studies 0, checked that in the thermody- 

namic limit both conditions agree as concern the zero-temperature effective action 
Eq. dO). 

The plan of the paper is as follows. In Section 2 we consider the Abelian constant 
chromomagnetic background field for the pure gauge SU(3) lattice theory at zero 
temperature. Section 3 is devoted to the study of both SU(2) and SU(3) lattice 
theories at finite temperature in presence of the external background field. Finally 
our conclusions are drawn in Section 4. 



2. SU(3) IN A CONSTANT CHROMOMAGNETIC FIELD 

We are interested in the case of a constant Abelian chromomagnetic field which in 
the continuum reads: 



AT\x) = A''^\x)6a,3 , (^) = h,2XiH 



:2.i] 
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The lattice links corresponding to A'^^^lx) can be evaluated from: 



[/^ = P exp \ ig I dt A^{x + tfi) 



(2.2) 



where P is the path-ordering operator, and = the A"'s being the Gell-Mann 

matrices. From Eqs. (|2.1|) and ( p.2|) we get: 

Ur\x) = f/|"*(f) = f/r*(x) = 1 



X) 



e 2 



■ gHxi 
' 2 











(2.3) 



Our Schrodinger functional Eq. ( |1.10 ) is defined on a lattice with the hypertorus 
geometry, for it is natural to impose that: 



U2{Xi, X2, Xs, Xi) = U2{Xi + Li, X2, X3, X^) 



(2.4) 



where Li is the lattice extension in the Xi direction (in lattice units). As a conse- 
quence the magnetic field H turns out to be quantized: 



2 



27r 



-n, 



ext ) 



(2.5) 



with next integer. 

According to the discussion of the previous Section, in evaluating the lattice func- 
tional integral Eq. ( p..lO| ) we impose that links belonging to the time slice X4 = are 
frozen to the configuration Eq. ( |2.3| ). Moreover we impose also that the links at the 
spatial boundaries are fixed according to Eq. ( p.3|) . In the continuum this last con- 
dition amounts to the usual requirement that the fluctuations over the background 
fields vanish at infinity. 

An alternative possibility is given by constraining the links belonging to the the 
time slice x^ = and those at the spatial boundaries to the condition 



f/fc(x) = f/r(f) 



1,2,3 



(2.6) 



while the links U^^x) are unconstrained. Note that with the condition Eq. ( |2.(j| ) the 
time-like plaquettes nearest the frozen hypersurface X4 = behave symmetrically in 
the update procedure. Moreover, in this way the Schrodinger functional Eq. ( |1.1(J| ) 
is the zero-temperature limit of the thermal partition functional Eq. ( |1.23|) . In our 
previous studies we checked that in the thermodynamic limit both conditions 
agree as the effective action is concerned. 

Our numerical simulations at zero temperature have been done on a lattice of size 
L1L2L3L4 with Li = L4 = 32, while the transverse size Li = L2 = L± has been 
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varied from Lj_ = 8 up to L± = 32. To avoid the problem of computing a partition 
function which is the exponential of an extensive quantity we consider the derivative 
of the density of effective action Eq. ( |1.12| ) with respect to P by taking next 

(i.e. gH) fixed. Indeed, we have 

e'[A 




dZ[0]' 



Z[0] dp 

-(^ E ^^^Trf/^-(^)) ' (2.7) 

\ X,^l>U I ^cxt 

where the subscripts on the average indicate the value of the external links at the 
boundaries, and the f/^i.(a;)'s are the plaquettes in the (/i, z/) plane. Actually, the 
contributions to £:'[y4^^*] due to the frozen time-slice at 0:4 = and the fixed links at 
the spatial boundaries must be subtracted. Accordingly, we define the derivative of 
the internal energy density: 



\ x'gA,/t>A / 



where A is the ensemble of the internal lattice sites which occupy the volume ^\nt- 
To implement the constraint at the boundaries in the numerical simulations we 
update only the internal links, i.e. the links U^{x) with x G A. We use the over- 
relaxed heat-bath algorithm to update the gauge configurations. Simulations have 
been performed by means of the APEIOO computer. Since we are measuring a local 
quantity such as the plaquette, a low statistics (from 1000 up to 5000 configurations) 
is required in order to get a good estimate of e[^^. 

In Figure 1 we display the derivative of the energy density normalized to the 
derivative of the external energy density: 

4xt = ^ [1 - cos(^)] = ^ [1 - cos(^ne.O] (2.9) 

versus P for Li = L4 = 32 and 8 < Lj_ < 32. From Figure 1 we see that, as in the 
SU(2) gauge theory both the peak and the perturbative plateau of e[^^ decrease 
by increasing L^. In order to perform the thermodynamic limit we introduce the 
scaling variable: 

x = ^, (2.10) 
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where 



is the magnetic length, and 



«h = a/|| [2.11] 



L.^ = ^it (2.12) 



is the lattice effective linear size. As in the SU(2) case ^ we try the scaling law: 

^-a g;nt(/^,^ext,I^cff) ^ ^^^^ _ ^2.13) 
^ext 

Indeed, from Figure 2 we see that our numerical data can be arranged on the scaling 
curve It is remarkable that the value of the exponent a = 1.5 in Eq. (|2.13| ) 



agrees with the one we found for the SU(2) gauge theory From Eq. ( |2.13| ) we 



can determine the infinite volume limit of the vacuum energy density eint- We have: 
lim £i„t(/?,next,i:cfr) =4xt / d(3' k{(3') lim ( =0 (2.14) 

in the whole range of (3. As a consequence, in the continuum limit (Leg — ^ oo,/5 — > 
oo) the SU(3) vacuum screens the external chromomagnetic Abelian field in accor- 
dance with the dual superconductivity scenario |l^ . 

3. FINITE TEMPERATURE BACKGROUND FIELD EFFECTIVE 
ACTION 

Let us consider the gauge systems in an external chromomagnetic Abelian field at 
finite temperature. According to the discussion in Section 1, we are interested in 

the thermal partition function Zt A'^^^ , Eq. ( |1.23| ). On the lattice the physical 



temperature Tphysicai is given by: 

T^T^ = Lt (3.1) 

physical 

where Lf = L^^ is the linear extension in the time direction. In order to approximate 
the thermodynamic limit, the spatial extension Lg should respect the relation: 

Ls:> Lt. (3.2) 

To this end we perform our numerical simulation on 32^ x Lt and 64^ x Lt lattices 
by imposing: 

^ < 4. (3.3) 
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At finite temperature the effective action is defined tlirougli tlie free energy: 



1 Zt /l" 

Obviously, in tfie case of constant external cliromagnetic field tlie relevant quantity 
is the density of effective action: 

^[^ext] ^ (3.5) 

As is well known, by increasing the temperature the pure gauge system undergoes 
the deconfinement phase transition. In the case of pure SU(N) gauge theories it is 
known that the expectation value of the Polyakov loop in the time direction 

^=^E^Trnf/4(f,a;4) (3.6) 

^ X 2:4=1 

is the order parameter for the deconfinement phase transition. As a preliminary step 
we look at the behavior of the temporal Polyakov loop versus the external applied 
field. We start with the SU(2) gauge system at /5 = 2.5 on 32^ x 5 lattice at zero ap- 
plied external field (i.e. next = 0) that is known to be in the deconfined phase of finite 
temperature SU(2). If the external field strength is increased the expectation value 
of the Polyakov loop is driven towards the value at zero temperature (see Fig. 3). 
It is worthwhile to stress that this last result is consistent with the dual supercon- 
ductor mechanism of confinement. Similar behavior has been reported by Authors 
of Refs. [lTI|,|r2| within a different approach. If we now consider the SU(2) gauge 



system at zero temperature in a constant Abelian chromomagnetic background field 
of fixed strength (next = 1) and increase the temperature, we find that the pertur- 
bative tail of the /^-derivative of the free energy density /4t(/3, next)/£ext increases 
with 1/Lt and tends towards the "classical" value //nt(/3, ''^ext)/£ext — 1 (^^^ Fig. 4). 
Therefore we may conclude that, as the temperature increases, there is no screening 
effect in the free energy density, confirming that the zero-temperature screening of 
the external field is related to the confinement. 

The knowledge of /4it(/9, 'n.ext)/£^ext finite temperature can be used to estimate 
the deconfinement temperature T^.. In Figure 5 we magnify the peak region for 
different values of Lt. We see clearly that the pseudocritical coupling [3*{Lt) depends 
on Lt. To determine the pseudocritical couplings we parametrize f[^t{P, Lt) near the 
peak as: 

fLtiP,Lt) ^ aijLt) 

4xt a2{Lt)[(3-f3*{LtW + l- 

We restrict the region near P*{Lt) until the fits Eq. ( |3.7|) give a reduced x^/d.o./. ~ 
1. 
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Having determined (3*{Lt) we estimate the deconfinement temperature as: 



_ 1 1 



where 



where bo = (llA^)/(487r^) and is the color number. In Figure 6 we display 



Tc/Aiatt for different temperatures. Following Ref. we linearly extrapolate to 
the continuum our data for T^/Aiatt- In this way we obtain the following estimate 
of the critical temperature in the continuum limit 

T 

28.36 ± 1.38. (3.10) 



A 



latt 



Equation ( p. 101) is to be compared with the continuum limit of the critical temper- 



ature available in the literature 



A 



T 

' = 24.38 ± 2.18. (3.11) 

latt 



Our result Eq. ( p.lQ ) agrees, within two standard deviations, with the result given 
in Eq. (|3.11|) . Note, however, that the small discrepancy could be a true dynami- 
cal effect due to the external chromomagnetic field. Indeed, the behaviour of the 
Polyakov loop versus the external background field displayed in Fig. 3 seems to sug- 
gest that the critical temperature does depend on the applied background field. For 
dimensional reasons one expects that: 

H. (3.12) 

To check the expected behaviour Eq. ( |3.12| ) we need to vary the external chromo- 
magnetic field. We plan to do this study in a future work. 

For the time being, let us turn to the SU(3) pure gauge lattice theory at finite 
temperature. In this paper we limit ourselves to present our determination of the 
critical deconfinement temperature. In Figure 7 we display the peak region of the 
internal free energy density for different values of Lt, together with the fits Eq. ( p.7|) . 
Having determined the pseudocritical couphngs, the deconfinement temperature can 
be obtained from Eq. (|3.^ ). Performing the linear extrapolation to the continuum 
limit we get: 

T 

= 20.86 ± 3.02. (3.13) 

Alatt 
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Our estimate Eq. (|3.13| ) is in fair agreement with the continuum hmit of the SU(3) 
critical temperature available in the literature [|T^] 




29.67 ± 5.47. 



(3.14) 



However, we note once more that our data do not exclude a dependence of the 
critical temperature on the external chromomagnetic field. 

4. CONCLUSIONS 

We have studied the non-perturbative dynamics of the vacuum of SU(2) and SU(3) 
lattice gauge theories by means of a gauge-invariant effective action defined using 
the lattice Schrodinger functional. 

At zero temperature our numerical results indicate that even for the more in- 
teresting case of the SU(3) theory, in the continuum limit L^s ^ oo, /9 ^ oo we 
have: 



so that the SU(3) vacuum screens completely the external chromomagnetic Abelian 
field. In other words, the continuum vacuum behaves as an Abelian magnetic con- 
densate medium in accordance with the dual superconductivity scenario. 

The intimate connection between the screening of the external background field 
and the confinement is corroborated by the finite temperature results. Indeed our 
numerical data show that the zero-temperature screening of the external field is 
removed by increasing the temperature. Moreover, at finite temperature it seems 
that confinement is restored by increasing the strength of the external applied field. 

At finite temperature we find that the /3-derivative of the free energy density 
behaves like a specific heat. From the peak position of the /^-derivative of the free 
energy density we obtained an estimate of the critical temperature Tc/Aiatt that 
extrapolates in the continuum limit to values in fair agreement with previous de- 
terminations in the literature. Moreover, our data are suggestive of a non-trivial 
dependence of the deconfinement critical temperature on the applied external chro- 
magnetic field. We deserve to a future work the investigation of this interesting 
possibility. 
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Figure 1: The /^-derivative of the internal energy density Eq. ( p.8| ) versus (3 for SU(3) 
on a L X x L lattice. Plotted data refer to L = 32 and different values of the 
transverse lattice size L i . 
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Figure 2: The scaling curve for SU(3) obtained by re-scaling all lattice data for 
e[^^{P,ne^t,Les) according to Eq. ( |2.13|) . 
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Figure 3: The absolute value of the Polyakov loop versus riext, for SU(2), on a 32^ x 5 
lattice at P — 2.5. 



14 



1.5 - 



X 

CO 



0.5 - 







"1 1 r 



"1 r 



"1 r 



"1 1 r 








o 




= 32 




□ 




= 8 




o 




= 7 




A 




= 6 




< 




= 5 





^< <<}<<}<<} 
^ A aaaaaaaaaaaaaaa 

□ ^ o o o 



oooooooooooooooo 



J I I I I I I I L 



J I I L 



4 

p 



8 



Figure 4: The /9-derivative of the free energy density Eq. (|3.5| ) for SU(2) on x L 
lattices, with Lg = 32. 
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Figure 5: The peak region for the /3-derivative of the free energy density Eq. (|3.5|) 
for SU(2). The data are taken on x Lf lattices, with Lg = 64. Sohd hnes are the 
best fits to the data using Eq. ( p.7|) . 
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Figure 6: Our SU(2) lattice data for Tc/Aiatt versus the temperature (circles). The 
full square is the continuum extrapolation of Ref. [l^ (see Eq. ( p.ll|) ). The dotted 
line is a linear fit to our data. 
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Figure 7: The peak region for the /3-derivative of the free energy density Eq. ( |3.5| ) 
for SU(3). The data are taken on x Lf lattices, with Lg = 64. Sohd hnes are the 
best fits to the data using Eq. ( p.7|) . 
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Figure 8: Our SU(3) lattice data for Tc/Aiatt versus the temperature (circles). The 
full square is the continuum extrapolation of Ref. [1^ (see Eq. ( p.l4| )). The dotted 
line is a linear fit to our data. 
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